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O " ABSTRACT 

We study the dependence of the galaxy size evolution on morphology, stellar mass 
and large scale environment for a sample of 298 group and 384 field quiescent early- 
type galaxies from the COSMOS survey, selected from z ~ 1 to the present, and with 
masses log(M/M Q ) > 10.5. 

J> | From a detailed morphological analysis we infer that ~ 80% of passive galaxies 

C*~) . with mass log{M /M©) > 10.5 have an early-type morphology and that this fraction 

0^ 1 does not evolve over the last 6 Gyr. However the relative abundance of lenticular and 

elliptical galaxies depends on stellar mass. Elliptical galaxies dominate only at the 
very high mass end - log(M /M©) > 11 - while SO galaxies dominate at lower stellar 
C-x ; masses - 10.5 < log{M/M Q ) < 11. 

<^> The galaxy size growth depends on galaxy mass and early-type galaxy morphol- 

ep ■ ogy, e.g., elliptical galaxies evolve differently than lenticular galaxies. At the low mass 

end - 10.5 < Log(M /M©) < 11, ellipticals do not show strong size growth from z ~ 1 
to the present (10% to 30% depending on the morphological classification). On the 
other end, massive ellipticals - log{M /M®) > 11.2 - approximately doubled their size. 
Interestingly, lenticular galaxies display different behavior: they appear more compact 
^ ■ on average and they do show a size growth of ~ 60% since z = 1 independent of stellar 

mass. We compare our results with state-of-the art semi-analytic models. 

While major and minor mergers can account for most of the galaxy size growth, 
we find that with present data and the theoretical uncertainties in the modeling we 
cannot state clear evidence favoring either merger or mass loss via quasar and/or 
stellar winds as the primary mechanism driving the evolution. 

The galaxy mass-size relation and the size growth do not depend on environment 
in the halo mass range explored in this work (field to group mass log(Mh/M & ) < 
14), i.e., group and field galaxies follow the same trends, which is at variance with 
predictions from current hierarchical models that show a clear dependence of size 
growth on halo mass for massive ellipticals - Zog(M*/M©) > 11.2. 

BCG and satellite galaxies lie on the same mass-size relation, at variance with 
predictions from hierarchical models, which predict that BCGs should have larger 
sizes than satellites because they experience more mergers in groups over the halo 
mass range probed. 

Key words: galaxies: evolution, galaxies: elliptical and lenticular, cD, galaxies: 
groups 
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1 INTRODUCTION 

The fact that massive quiescent galaxies experienced a 
strong size evolution in the last 10 Gyrs is now a com- 
monly accepted picture since first works on thi s topi c 
were published ((Daddi et all 120051 ; iTruiillo et ail l2006h . 
Many independent groups using different datasets and 
selections have come up with similar conclusions, i.e. 
massive galaxies roughly doubled their size from z ~ 

1 and probably multiplied it by 3 ~ 5 from z ~ 

2 (e.g. Ivan der Wei etafl 120081; Ivan Dokkum et all |200S|; 



Buitrago et all 120081 ; iDamianov et all l201ll ; ICimatti et al.1 
20121 ) even though there might be a populati on with larger 



sizes already in place at high redshift (e.g iMancini et al.l 
l20ld : | aracco et al.l 1201 if ) As several works pointed out 
since the first publications came out, the result could be 
biased by a wrong estimate of stellar masses (which i s usu - 
ally done through SED fitting, e.g. IRaichoor et al.1 1201 ll ) 
and/or an under-estimate of galaxy sizes at very high red- 
shifts since surface brightness dimmin g could cause the loss 
of the very outer par ts of galaxies (e.g. lBezanson et al.ll2009l ; 
IHopkins et all 12009). Independent measurements based on 
dynamical masses have howe ver confirmed the compact- 
ness of several objects ( e.g. [Martinez-Manso et al.l 1201 ll : 
van de Sande et~alTl201ll ; [n cwman et al.1 12012| )" and it now 
seems clear that there is indeed a population of compact 
objects at high redshift. Some works also pointed out that 
we cannot exclude that the current galaxy selections might 
be biased and we are missing those compact objects in the 
nearby univer se, or selecting the most compact objects at 
high redshift (|Valentinuzzi et all 2010l) . This would be at 
variance with the work of ITruiillo et all (|2009l ) that showed 
that there is not a strong evidence for a a significant fraction 
of compact objects in the local Universe. 

From the theoretical point of view, two ma in mecha- 
nisms have been proposed to puff up galaxies. IFan et al.l 
(|2008l ) proposed AGN and supernovae feedback as the main 
resp onsible of galaxy e xpansion while IHopkins" et al.l (|2006l ) 
and INaab et al.1 (|2009l ) argued that minor dry mergers are 
the most efficient way to grow sizes since they affect the 
outer parts of the galaxy without significantly modifying the 
stella r mass nor the star formation (see also IShankar et al.l 

Hon]). 

Several recent observational works have tried to dis- 
entangle the two scenarios. As a result, minor dry merg- 
ing seems to emerge as the most plausible explanation for 
the size evolutio n at least from z ~ 1 even though there is 
still some debate. [Newman et al. (2012) showed by carefully 
counting small com panions in very de ep images from the 
CANDELS survey (jGrogin et al.l l201ll ) that minor merg- 
ers are roughly enough to account for the size evolution 
from z ~ 1 (provided that a short merger timescale is 
assumed) and a combination of minor merging and star 
formation quenching can explain the galaxy growth from 
z ~ 2. Using morphological merging indicators, Bluck et al.l 
(|201ll ) reached a similar conclusion and ev en argue that 
the p roblem is close to be solved (see also IMcLure et al.l 
||2012|) for similar consideration s). On the other hand, 
ILopez-Saniuan et al.1 (|201ll . l2012h point out that for galax- 
ies with masses log(M /Mq) > 11 minor mergers are not 
the only process responsible for size growth since z ~ 1. 
These authors propose a scenario for which minor and ma- 



jor mergers contribute to ~ 55%, while the remaining ~ 45% 
to ~ 25% is due to other processes and specially to younger 
galaxies (hence larger) arriving at later epochs (progenitor 
bias). 

When measuring the galaxy mass-size relation and 
mass growth, different works use different sample selections 
though and this might lead to different conclusions. It is very 
rare in fact to find two works dealing with size evolution 
which apply the same criteria to select their galaxy sam- 
ple and the same methodology to a nalyze the data . Some 
are based on star-for mation only ([Papovich et al. 2012T1 
others on morphology (|Raichoor et al.ll2012l; Cooper et al 



200S 



2011b) or on a combination of both (e.g. Ivan der Wei et al 



Damianov et al.ll201ll ). Finally, many works combine 



a stellar mass selection with a quiescence criterion (e.g. 
Ivan Dokkum et al.|[2008l ; lNewman et al.|[2012l) and generally 
the m ass c uts are not always the same. See lDamianov et al.l 
(|201ll ) and ICimatti et"aT1 (|2012l ) for two different compila- 
tions of recent results. 

Another recent point of discussion in the literature 
is the role of environment . Recent observationa l studies 
show controversial results. IRaichoor et~al] (|2012l ) studied 
the mass size relation for morphologically selected early- 
type galaxies at z ~ 1.2 in three different environments 
(field, cluster, groups) and find that, on average, for masses 
10 < log(M /Mq) < 11.5 cluster galaxies appear to be 
smaller at fixed stellar mass than field galaxies. Interest- 
ingly, in the same ste llar mass range but at lower redshift, 
ICooper et all ll2011bl) find ex actly the opposite trend using 
DEEP3 (|Cooper et al.ll2011al ). Larger sizes in the cluster en- 
vironment are also observed at z = 1.62 bv lPapovich et all 
|2012l ) with CANDELS data for passive galaxies with stel- 
lar masses larger than log(M /Mq) ~ 10.5. Howeve r, other 
two works (|Maltbv et alj ' l20ld : iRettura et~alll20ld ) do not 
find any trend with environment at z < 0.4 and z ~ 1.2 re- 
spectively. The differences between these works are somehow 
puzzling but might come from the different sample selections 
(e.g.: based on color or Sersic index vs visual morphology 
classification) and/or the way environment is measured (lo- 
cal vs. global) and /or low statistics (high redshift results 
indeed rely on some tens of galaxies). 

We will address here these questions by carefully se- 
lecting a statistically significant sample of galaxies from the 
recently published COSMOS X-ray detected grou p galaxy 
sample |George et al.|[201ll ; iFinoguenov et al.|[2007l ) as com- 
pared to the field. Within this sample we will dissect the 
properties of passive galaxies and look at the effects on the 
size growth of the galaxies with different morphology in dif- 
ferent environments. 

The paper is organized as follows: in section [2] we de- 
scribe the datasets used in this work and the derived quanti- 
ties (morphologies, sizes and stellar masses) and in sections[3] 
and U we present and discuss our main results. Through- 
out the paper we consider a standard ACDM cosmology 
(fiaf = 0.3, n A = 0.7). 



2 DATA AND ANALYSIS 
2.1 Datasets 

We use two samples of galaxies from the COSMOS survey 
belonging to two different environments: groups and field. 
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The group sample is composed of groups in the COS- 
MOS fie2d_tha^Jiave_tee^i_detected as extended X-ray emit- 
ters ()Finoguenov et all 120071 ) , which is a strong signature 
of virialized structures, and have several spectroscopic con- 
firmed members. The full sample contains groups with halo 
masses from M2 00C /Mp> ~ 10 13 to ~ 10 14 as measured 



by weak-lensing ( Leauthaud et al]|201Ch and spans the red- 
shift range 0.2 < z < 1.0. Group members have been 
selected based on photometric redshifts deriv ed from the 
exten sive COSMOS multi-wavelength imaging (|llbert et al.l 
I2009TI. More det a ils abo ut the galaxy selection can be found 



George et al.l (|201lh who carried out a careful analysis 



of potential biases and contaminations and found that the 
galaxy member sample has a purity of 84% and a complete- 
ness of 92% within 0.57?2ooc- For this work, we select all 
galaxies within -R200C with a probability of being in the 
group greater than 0.5 and Isu(AB) < 24 mag - this last 
cut required for size estimates as explained in section 12.21 
This selection results in a group sample of 3146 galaxies dis- 
tributed over 70 groups. We expect th en a contamination of 
» 30%, and a completeness of « 90% (jGeorge et al.ll201ll ). 

Field galaxies are selected in the COSMOS field in the 
same redshift range as group members and with the same 
magnitude cut but making sure that they do not belong 
to any detected gro up (e.g. with GROUPid = — 1 in the 
iGeorge et al.1 <|201lh group catalog). We assume that these 
galaxies lie in a dark matter halo of mass M200C / 'Mq < 
10 13 ' 2 , otherwise they shoul d have been detecte d as groups 
members - see Fig 1 from IGeorge et al.1 (|201ll ). The field 
sample contains 3760 galaxies. 



2.2 Sizes and masses 

Sizes of all galaxies have been estimated on the COSMOS 
HST/ ACS F814W image s (Mosaic v2.0 iKoekemoer et al.1 
l2007h using GALAPAGOS l|Barden et all l2012h which is an 
IDL based pipeline to run Sextrac tor |Bertin fc Arnoutsl 
119961 ') and Galfit (|Peng et al.ll2002r i together. We basically 
fit every galaxy in the field and in the group sample with 
a 2D Sersic profile <|Sersidll96^ ) using the d efault GALAPA- 
GOS p arameters as described and tested in lHaussler et al.l 
l|2007h . Our number of failures (i.e. fits that do not con- 



verge) is less than 5%. The poi nt spread function ( PSF) 
used for the fitting is taken from iRhodes et alj (|2007l ) who 
computed spatially-varying model PSFs for the COSMOS 
survey taking into account variations in the effective HST 
focus positions. For this work we used a single PSF (esti- 
mated at the average focus position for the COSMOS survey 
(A = — 2/im)) for all the galaxies. 

The reliability of our size estimates is estimated by 
placing mock galaxies (1 < n < 8, 0.1 < r e < 1.5 , 
17 < 7 < 24 mag) in a real background. In our simulations 
we also explore the possible biases due to local over densities 
by dropping the mock galaxies in the same positions as in 
real high redshift clusters. We find that our size measure- 
ments are unbiased (| < (r e , out — r et i n )/r e:in > | < 0.1) with 
a reasonable scatter ( ^JVar[{r e , out — r e ,in)/r e ,in\ < 0.2) up 
to hn{AB) < 24 < 1.0 arcsec and /xsiw < 24 

mag x arcsec -2 (Fig. [T] and Delaye et aL, in prep). All 
galaxies in our sample have surface brightness brighter than 



M814W = 24 mag x arcsec~ so our size estimates are reliable 
according to the simulations. 

Since we are using a unique wavelength to estimate 
sizes up to z ~ 1, the rest-frame filter of our size estimates 
moves from the r band at z ~ 0.2 to the U band at z ~ 1 
which could create an artificial difference when considering 
the evolution between low and high redshift data. However, 
as shown in section [3] and in figure [7] the size evolution we 
measure in our sample for early-type galaxies is fully consis- 
te nt with previous publis hed results, specially with the ones 
of E ewman et al.l (|2012T l who computed s izes in a unique 
rest-fr ame filter. Moreove r , rece nt work bv lDamianov et al.l 
(|201ll ) and lCassata et al.l (|201ll) shows that sizes in the ul- 
traviolet and optica l rest-frame strongly correlate (see also 
ICimatti et al.l [2012h . Therefore we do not expect a signifi- 
cant bias due to this issue. 

Throughout the paper we will use circularized effec tive 
radii as primary size estimator, i.e. r" rc = r{ lt x y^b/a. 

Stellar masses are computed through spectral en- 
ergy distribu tion (SEP) f i tting using the Bayesian code 



described in iBundv et al.1 (|2006l ) (KB06) which empl 



oys 



iBruzual fc Charlotl (l2003h fBC03) models with a Chabrier 
IMF, and published in IGeorge et all J201l|) catalogs. More 
detail s can be found in section 2.3 of [Leauthaud et al.l 
l|20ld ). 



2.3 Morphology 

Galaxy morphologies in the group and field samples are also 
computed on the HST/ACS F814W images with GAlSVM, 
a code for automated morphological classification based on 
support ve ctor machines specially designed fo r high red- 
shift data (|Huertas-Companv et all 120081 . 12001 ). The code 
requires a visually classified training set which is used 
to simulate galaxies at higher redshift. The training set 
used i n this work is a combination of the lNair fc Abraham! 
|2010f) sample of ~ 14000 galaxies and the EFIGI project 
|Baillard et ail 120111 ) both from the Sloan Digital Sky Sur- 
vey DR4. Our final training set therefore contains ~ 20.000 
galaxiesQ] We follow the bay esian approach presented in 
iHuertas-Companv et al.l (|201 ih to associate to every galaxy 
4 probabilities of being i n 4 morphological classes as defined 
in the local universe bv lNair fc Abraham! (|20i0r i, i.e. ellip- 
ticals (E0, E+), lenticulars (SO—, SO, S0+, S0/a), early spi- 
rals (Sa, Sab, Sb, Sbc) or late spirals (Sc, Sd, Sdm, Sm, 7m): 
P(E),P(S0),P(Sab),P(Scd). Errors in probabilities are 
computed by bootstrapping, i.e. we repeat the classi- 
fication 10 times with randomly selected training sets 
from the main sample and keep the average proba- 
bility for each galaxy as the final classification (see 
IHuertas-Companv et al.l l201ll for more details). We then 
select as elliptical and lenticular galaxies the galaxies for 
which max(P(E),P(S0),P(Sab),P(Scd)) = P(E) and 
max(P(E),P(S0),P(Sab),P(Scd)) = P(S0) respectively. 
Early-type galaxies are then defined as the combination of 
both populations. With this selection, the completeness for 
early-type galaxies is ~ 95% and the contamination rate is 



1 An updated and stable version of the code as well 
as the training set used for this work arc available at 
http://gepicom04.obspm.fr/galSVM/Home.html 
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Figure 1. Results of simulations to assess the accuracy of our size estimates. Left panels show the difference between the input and 
output magnitudes (top left), effective radii (middle left) and Sersic index (bottom left) for individual simulated galaxies as a function 
of surface brightness. Right panels show the average bias and dispersions for different surface brightness bins. 



expected to be less than 7% as stated by detailed compar- 
isons with visual morpholog ical classifications at z ~ 1.3 
from (|Mei et al. 1 120091 . |2012| ) in the Lynx super cluster at 
2 = 1.26. Our classification is therefore very close to a vi- 
sual classification. This point is also confirmed by the fact 
that the axis ratio distribution of our early-type sam ple is 
very close to the one reported bv lBuitrago et all |201ll ) from 
a visually classified sample (see Fig. [2] and their figure 2). 

To compare our morphological classification with those 
often used in the literature, we estimated the contamination 
that would suffer a selection of early-type galaxies, based 
on a simple n > 2.5 cut. We find that such a sample would 
be c ontaminated by approximately 50% of early spirals (see 
also iMei et al.l [20121 ). We refer to section [3] for a detailed 
analysis of the implications of such a selection. 

Separating ellipticals from lenticulars is extremely chal- 
lenging even by eye, so our classification is necessarily more 
contaminated. Simulations show that we are close to 30% 
uncertainties, similar to those observed in o ther works (e.g. 
by v isual morphological classifications by IPostman et al.l 
120051 ). Our visually trained probabilistic approach for galaxy 
classification allows us to distinguish galaxies based on a 
quantitative separation in the parameter space that corre- 
sponds to that usually used to separate SOs and ellipticals in 
the local Universe (e.g. by visual classification). Indeed, the 
stellar mass and axis ratios distributions of the two popula- 
tions shown in Fig. [2] present different behaviors, as stated 
by Kuiper tests (P < 0.6), confirming that we are seeing 
separate populations and that the classification is not ran- 
dom. Ellipticals appear rounder, more massive on average, 
with slightly higher Sersic indices and larger sizes. In order 
to double check our separation between lenticulars and ellip- 



ticals we performed two ad ditional tests. First , we t ook the 
visually classified sample of iNair fc Abraham] (|2010h in the 
SDSS, cross-correla ted it with the Sersic decompositions re- 
cently published by Si mard et al.l (|201ll ) and looked at the 
same properties we studied for our high redshift sample. Re- 
sults are shown in Fig. We do clearly observe the same 
trends, i.e. lenticulars are also more elongated and slightly 
more compact than ellipticals, even though it is less pro- 
nounced than in the high redshift sample. This test confirms 
that our E/S0 classification is consistent with the local vi- 
sual classification, as expected. As a second step, two of us 
(MHC and SM) made visual classifications of our two auto- 
matically defined classes. We find that the two classification 
of E/SOs agree at ~ 80% at z < 0.5 and ~ 70% at z > 0.5 
which is fully consistent with the results from simulations 
and also with the di fferences found betwe en independent hu- 
man classifiers (e.g. lPostman et al.ll2005l ). However, in order 
to make sure that our results are not biased because of the 
automated method used, we double check our main results 
using visual classifications (see appendix O and comment 
in our analysis in subsequent sections whenever there is a 
significant difference. 



We deduce that the population of elliptical galaxies se- 
lected through our method is dominated by pure bulge sys- 
tems whereas an important fraction of what we call lentic- 
ulars have a disk component (but not observed arms) even 
if they are still bulge dominated. We also notice that a sim- 
ple Sersic based selection of any kind would not allow the 
separation between these two populations. 
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Figure 2. Stellar mass (top left panel), axis ratio (top right panel), Sersic index (bottom left panel) and size distributions (bottom 
right panel) of morphologically classified ellipticals (solid line) and SOs (dashed line) with the automated procedure described in the text. 
Error bars are Poissonian. SOs and ellipticals have different parameter distributions, which proves that the algorithm indeed separates 
two different populations of galaxies and does it in a quantitative and reproducible way. 



3 RESULTS 



3.1 Sample selection 



We investigate in this section how the sample selection af- 
fects the observed mass-size relation and size evolution of 
selected galaxies. We compare first three selections usually 
found in other works: 

(i) a passive selection. This selection includes all quiescent 
galaxies with log(M /Mq) > 10.5. Throughout this work, 
we will select quiescent or passive galaxies based on the 
M (NUV)-M(R) dust c orrected rest-frame color as computed 
bvlllbert et all (|2009l ) by means of SED fitting. Illbert et af] 
(|2009h have shown that a M(NUV) - M(R) > 3.5 selec- 
tion results in a g ood separation of pas sive and star-forming 
galaxies. See also iGiodini et al.l (|2012f ) for a discussion on 



this selection on the same sample as the one used for this 
work. 

(ii) a Sersic based selection, i.e. we apply a simple n > 2.5 
cut as usually done as well as the same stellar mass selection 
(log(M/M Q ) > 10.5). 

(iii) an ETG selection in which we select all early type 
galaxies (ellipticals and lenticulars) with log(M/M®) > 10.5 
from our morphological classification detailed in section [2~3l 
independently of the star formation activity. 

We show in Figs. [4] to [6] the mass size relations for sam- 
ples (i), (ii) and (iii) respectively in different environments 
and redshifts and we summarize the best fit parameters com- 
puted through a standard chi square minimization in table[T] 
Also are shown in the figures, with different symbols the 
morphologies of the objects belonging to the given selection 
as well as the star formation activity (quiescent or active) 
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Figure 3. Same plot as figure [2] but for a local sample from the SDSS (Simard et al. 2011). The same trends as for the high rcdshift 
sample are observed. 



color coded for selection (ii) and (iii) . We also show the po- 
sition of the brightest group galaxies (BGGs), defined as 
the most massive galaxies wi thin an NFW scale radius of 
the X-ray position, from the iGeorge et al.l f|201 lh catalog. 
Notice that some of them are morphologically classified as 
early-spirals. This is probably because of the extended halo 
usually found in these objects which could be interpreted as 
a disk component by the automated algorithms. 

Despite the fact that we see some differences between 
the different selections, the fit parameters are roughly con- 
sistent within la. Since a Sersic index based selection con- 
tains an important fraction of star-forming galaxies (~ 50%) 
(a population which presents, on average, higher sizes than 
ETGs; e.g. Mei et al. 2012), this population tends to increase 
the scatter of the relation. 

A similar behavior is observed if we select galaxies just 
based on the morphology (ETGs, fig. O ■ The scatter of the 
relation might be increased by the presence of a star-forming 
population of galaxies with early-type morphology. Even if 



the number densities of these object s are small in the local 
universe (e.g 
at high z (e.g 



Kannappan et al. 20091) they tend to increase 



Huertas-Companv et alj|2010l ) and might have 



a consequence on the size measurements. 

We quantify the possible effects of galaxy selection in 
Fig. [7] where we plot the mass-normalized radius (R e x 
(10 11 Mq/M*)) evolution for different galaxy selections 
(ETG, passive ETGs, n > 2.5, passive galaxies) with 
log{M /Mq) > 10.7 (the mass selection for this fig ure is 
chosen to match the one from [Newman et alj 12012 ). We 
observe that surprisingly, the differences between the dif- 
ferent selection in terms of size evoluti on are negligible and 
also consistent with recent results (i.e . Cimatti et al.l [20121 ; 
iDamianov et al.ll201ll ; iNewman etal] |2012h . We can con- 
clude that these different selections usually found in the lit- 
erature lead to similar results. 

In this work we have an additional ingredient though, 
usually lacking in previous published results which is the 
morphological dissection of passive galaxies as explained 
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Figure 4. Mass-size relation of passive galaxies in groups (left column) and field environments (right column). Filled circles are ellipticals, 
empty squares are lenticulars and stars are late- type galaxies (including early and late- type spirals). Diamonds show the position of BGGs. 
We also show on the bottom right corner the typical error bars for sizes and stellar masses. The solid line shows the local relation as 
measured by Bernardi et al. (2010) and dashcd-dotted lines are the best fits to the whole passive population. The distribution in the 
mass-size plane of the different morphological types are not the same. A pure star-formation selection does not ensure that we are 
studying just bulge growth but also galaxies with a disk component which might follow a different evolutionary path. 
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Figure 5. Mass-size relation of n > 2.5 galaxies in groups (left column) and field environments (right column). Filled circles are 
ellipticals, empty squares are lenticulars and stars are late-type galaxies (including early and late-type spirals). Diamonds show the 
position of BGGs. Blue/red symbols show active/passive galaxies respectively. We also show on the bottom right corner the typical error 
bars for sizes and stellar masses. The solid line shows the local relation as measured by Bernardi et al. (2010) and dashed-dotted lines 
are the best fits to the whole n > 2.5 population. A Sersic index based selection contains an important fraction of star-forming galaxies 
which tend to increase the scatter of the relation. The effect seems to be more pronounced in the field population so it could yield to 
observed differences in the evolution of the mass-size relation with environment 
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Figure 6. Mass-size relation of early-type galaxies in groups (left column) and field environments (right column). Filled circles are 
ellipticals and empty squares arc lcnticulars. Diamonds show the position of BGGs. Blue/red symbols show active/passive galaxies 
respectively. We also show on the bottom right corner the typical error bars for sizes and stellar masses. The solid line shows the local 
relation as measured by Bernardi et al. (2010) and dashed-dotted lines are the best fits to the whole early-type population. 
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Figure 8. Example stamps of passive galaxies of 4 morphological classes. First row: ellipticals, second row: SOs, third row: early-spirals 
and last row: late-spirals. Stamp sizes are: 1 X 1 . 
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Figure 9. Morphological fractions of passive galaxies as a function of stellar mass and redshift. Red circles: ellipticals, cyan triangles: 
SOs, green squares: early-spirals, blue diamonds: late-spirals. The black dashed line shows the fraction of early-type galaxies (ellipticals 
and SOs together). The uncertainties are calculated following Gchrels (1986; see Section 3 for binomial statistics; see also Mei et al. 
2009). These approximations apply even when ratios of different events are calculated from small numbers, and yield the lower and upper 
limits of a binomial distribution within the 84% confidence limit, corresponding to la. Note that using this conservative approach our 
uncertainties are slightly overestimated (Cameron et al. 2011). 



in section 12.31 As a matter of fact, the population of 
passive galaxies is not a homogeneous population of ob- 
jects. We show indeed some example stamps of massive 
(log(M / Mq) > 10.5) passive galaxies with different mor- 
phologies in Fig. [8] It is easy to notice by simple visual in- 
spection that not all passive galaxies are bulge dominated. 
The relative abundance of each morphological type as a 
function of stellar mass is quantified in Fig. [9] for all our 
passive galaxy sample. 

As expected, early-type morphologies dominate the 
population of passive galaxies up to z ~ 1 at all stellar 
masses, being about 80% of the total population at all red- 
shifts. The remaining 20% is populated by early-type spirals 
while late-type spiral fractions are negligible. At all redshifts 
up to z ~ 1, elliptical galaxies dominate the ETG popula- 



tion at masses log(M / Mq) > 11 — 11.2. At z < 0.5, galaxies 
with masses log{M / Mq) < 11 are around half ellipticals 
and half lenticular, while lenticulars dominate the low-mass 
fractions at z > 0.5. In other words, an important fraction of 
passive galaxies from z ~ 1 (if not the majority) could have a 
disk component (see alsolBundv et al 1 l20TOl . iMei et alJlioil 
van der Wei et al]|201ll ). We also show in appendix, a ver- 
sion of figure [9] obtained with the two independent visual 
classifications performed in this work. One of the classi- 
fiers (SM) finds visually more Sas than SOs with respect 
to the automated classification, specially at high redshift. 
However the separation between galaxies with disk (Sas and 
SOs) and without disk (ellipticals) is very similar to what 
is obtained with the automated classification, except that 
elliptical galaxies are not dominant at the high mass end in 
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Table 1. Power-law fitting parameters to the mass-size relation (r e = (-jnjQ ) a + b) for different galaxy selections and environments. Pall 
stands for Passive all (sample (i) from the text), PEs is Passive ellipticals, PSOs means Passive lenticulars, n > 2.5 are galaxies from 
a Sersic based selection (sample (ii)) and ETGs stands for early-type galaxies independently of star formation (sample (iii)). N is the 
number of objects in each bin with log(M* /Mq) > 10.5. 



the highest redshift bin, which is also found by the second 
classifier (MHC). 

Studying the size evolution of passive galaxies all to- 
gether (as often done in several works), mixes not only ETGs 
with early-spirals (Sa-b), but also ETG different morpholog- 
ical populations (e.g. ellipticals and lenticulars) for which 
the evolution is not necessarily of the same nature. 

Differences in the mass-size relation of ellipticals and 
lenticulars are visible in Fig. [4] to [6] and Table [1] SO galaxies 
tend to be not only less massive but also more compact with 
respect to ellipticals. 

In the next sections, we will focus on the sample of pas- 
sive galaxies (sample (i)) and will quantify these differences 
by studying the dependence of the size evolution of passive 
galaxies on different morphological types (namely ellipticals 
and lenticulars). The final sample contains 298 group and 
384 field galaxies. 

3.2 Size growth of passive galaxies with different 
morphologies from z ~ 1 

Throughout this sec tion we normal i ze siz es using the local 
reference derived bv lBernardi et al.1 (|2010l 'l on the SDSS us- 
ing a clean sample of elliptical galaxies (excluding SOs). We 
payed special attention in properly calibrating the local ref- 
erence to be as consistent as possible with our high redshift 



measurements. We checked for that purpose that our low 
redshift (z < 0.3) data are consistent at la with their lo- 
cal mass-size relation which confirms that both relations are 
well calibrated. Sizes in each redshift bin are computed by 
fitting a 1-D gaussian function to the size ratio distribu- 
tions and keeping the position of the peak as the adopted 
value. Results remain however unchanged if we use instead 
a classical median or a three sigma-clipped average. Uncer- 
tainties are computed through bootstrapping, i.e. we repeat 
the computation of each value 1000 times removing one ele- 
ment each time and compute the error as the I — a error of 
all the measurements. 

3.2.1 Ellipticals 

In fig. 1101 we show the size evolution of ellipticals as com- 
pared to lenticulars in two mass bins (10.5 < Log(M / Mq) < 
11 and 11.0 < Log(M/M Q ) < 11.5). The upper limit is cho- 
sen because there are almost no SOs with masses greater 
then 10 1 solar masses (see fig. [2} and we want to compare 
the evolution of the two morphological types in the same 
mass ranges. 

Ellipticals in the mass range 11.0 < Log{M / Mq) < 11.5 
experienced a ~ 40% ± 10% size growth from z ~ 1 to 
present. We find an a value of a = 0.8 ± 0.28 when fit- 
ting the evolution with a power-law (r e cx (1 + z)~ a ). How- 



The dependence of the ETG mass-size relation on mass, morphology and environment from zr~>l to present 



ever, the most relevant feature is that, as shown in the left 
panel of fig. 1101 below ~ 10 11 solar masses, elliptical galax- 
ies do not experience a significant size growth, e.g their size 
does not evolve in a significant way with respect to local 
ellipticals from the S PSS (r e oc (1 + z )-°- 34 ± 017 ) ( see also 
iRaichoor et al.l [20121 for similar results at higher redshift). 
The mass dependence is even more clear in figure [11] in 
which we show in the left column, the observed size evo- 
lution for ellipticals in three bins of increasing stellar mass 
(10.5 < log(M/M Q ) < 11, 11 < log(M/M Q ) < 11.2 and 
11.2 < log(M/M Q ) < 12). The derived a value parameter- 
izing the evolution is a = 0.34 ± 0.17, a = 0.63 ± 0.18 
and a = 0.98 ± 0.18 from low to high mass, showing a 
clear increasing trend with stellar mass. In appendix [X] we 
show that when using visual morphological classifications 
the evolution in the low mass bin is stronger (~ 30%) mak- 
ing the difference between the two first mass bins less pro- 
nounced. The trend is still preserved though (see table |A"Tt . 
specially between the two last bins. One possible source of 
error though is that the stellar mass bins used are compara- 
ble to the expected error of the stellar mass (i.e. 0.2 dex), so 
the different behaviors found can be affected by contamina- 
tions of objects with lower/higher stellar masses. We have 
run Monte Carlo simulations to check if the trends found are 
preserved with typical errors expected on the stellar mass. 
For that purpose, we added to every stellar mass a random 
shift within 3a of the expected error in stellar mass and 
recomputed the median sizes 1000 times. The values found 
are consistent at la level so we do not expect a significant 
contribution of this effect in our measurements. 

Interestingly, this mass dependence is less pronounced 
when studying all ETGs or passive galaxies as a whole (right 
column of fig. Ill|) . The best fit relations for this whole sam- 
ple are r e oc (1 + z )-i-0i±0.23 for ET q s witn ste llar masses 
10.5 < Log{M / Mq) < 11 and r e oc (1 + z )- 1 - 21 ± - 22 an d 
r e oc (1 + 2 )- 119±0 18 for 11 < Log(M/M Q ) < 11.2 and 
11.2 < Log(M /Mq) < 12 respectively which are fully con- 
sistent within la. 

The mass dependence has been discussed in the recent 
liter ature providing d ifferent results. IWilliams et al.l (|2010l ) 
and lRvan et alj (|2012T ) measured a mass dependence similar 
to the one reported here for ellipticals, i.e. with the evolution 
bei ng stronger at higher st ellar masses, while other works 
like iDamianov et al.l (|201lf) suggest the slope of the mass- 
size relation is mass independent. We show here that the re- 
sults might depend on how the selection is performed. When 
selecting pure passive bulges a mas s dependent ev o lution 
seems to emerge. As a matter of fact. I Williams et al.l (|2010l ) 
select passive galaxies based on the specific star formation 
rate (instead of the red sequence for other works). They 
might be removing from their sample more disky galaxies 
(also removed from our elliptical sample and not from an 
ETG sample) which could explain that we find similar re- 
sults. Notice also that for galaxies with log(M r /Mo) > 11.2, 
the contribution from disky galaxies is almost zero and 
therefore selecting ETGs or ellipticals has almost no effect in 
the size evolution. The strong size evolution found for these 
galaxies is therefore robust to galaxy selections and will be 
discussed in detail in section [4] 



3.2.2 Lenticulars and early spirals 

The bottom panels of Fig. [10] show the size evolution of 
lenticular galaxies in two mass bins normalized to the same 
local relation. They appear on average more compact than 
ellipticals, in a similar way as in the local Universe, as shown 
in sectionl2.3land confirm ing the fact already pointed out by 
Ivan der Wei et ail (|201lf ) that the most compact galaxies at 
high redshift have a disk component (this effect is preserved 
even when not circularizing the radii). Since we normalize 
with the same local relation as for ellipticals, we also plot, 
for reference, the ratio of lenticular sizes over those of the 
ellipticals in the SDSS (notice that all fractions are normal- 
ized to the local SDSS elliptical sizes). 

What is interesting is that the size growth of lentic- 
ulars does not seem to depend significantly on stellar mass 
given the large uncertainties, contrary to the behavior shown 
for the elliptical population: r e oc (1 + z )-°- 67±0 ' S0 an d 
r e oc (1 + 2;)- 102 ± 25 f or i ow m ass and massive lenticulars, 
respectively. 

In order to further investigate what is driving the 
size increase in these galaxies, we made a simple exer- 
cise which is to compare the observed evolution to the 
one expected in a star-forming disk dominated popula- 
tion. We selected for that purpose, galaxies in our sam- 
ple with M(NUV) - M(R) < 3.5 and a spiral-like mor- 
phology (max(P(E),P(S0),P(Sab),P(Scd)) = P{Sab) or 
max(P(E),P(S0),P(Sab),P(Scd)) = P(Scd)), computed 
the size growth as a function of redshift as we did for the 
SOs. Since the overal l sizes of late-t ype galaxies are larger 
than ETGs (see also lMei et ai1l2012r i, we divide the overall 
relation by a factor of 1.3 so that the values at z = agree 
(see fig. UOp. We just plot the relation for the low mass end 
because there are not enough high mass late-type to derive 
a reliable relation at high masses. 

Even though the normalization for disks is still slightly 
higher, the trends of the size growth for the two populations 
are very similar (r e oc (1 + z )-° 67 ± a3 ° f or lenticulars and 
r e oc (1 + z) _0 ' 54±0 ' 44 for star forming disks) which suggests 
that the size growth of a large part of the early-type galaxies 
and that of the late-type galaxies evolve at similar rates. 
If we add to the passive SO population, also passive early- 
spirals (max(P(E),P(S0),P(Sab),P(Scd)) = P{Sab) and 
M(NUV) - M(R) > 3.5) the trend is preserved. 

3.3 Environment 

We now study the effect of environment on the mass-size 
relation of galaxies paying special attention to the uncer- 
tainties due_to_th«_g_alaxy_SMnple selection. 

In IRaichoor etail i|2012h . in a sample dominated by 
galaxies with 10 < log(M / 'Mq) < 11, we have noticed that 
while differences in median/average sizes cannot distinguish 
early-type populations in different environments, a test on 
both the mean and scatter of their distribution can point 
out environmental differences. When comparing cluster and 
field early-type mass/size distribution, we found that field 
ETGs (defined visually as ellipticals and SOs) at z ~ 1.3 
show a mass-size relation similar to the local one (e.g. their 
mass-size relation did not evolve overall from z ~ 1.3 to 
present) and a larger scatter with respect to cluster ETGs, 
that are also on average more compact that cluster ETGs in 
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Figure 10. Size evolution of passive elliptical (top) and SO (bottom) galaxies in two mass bins. Left panel: 10.5 < log(M / Mq) < 11, 
right panel: log(M / Mq) > 11.2. Dottcd-dashed lines are the best fit lines and red line in the bottom left panel is the measured evolution 
of star- forming galaxies (see text for details). Numbers indicate the number of objects in each redshift bin. Error bars are the scatter of 
the distributions. SDSS points are computed using the catalogs from Nair & Abarahm (2010) and Simard et al. (2011). 
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Figure 11. Size evolution of elliptical (left column) and early-type passive galaxies (right column) in three mass bins. Top row: 10.5 < 
log(M/M Q ) < 11, middle row: 11.0 < log(M/M Q ) > 11.2 and bottom row: 11.2 < log(M/M Q ) > 12. Numbers indicate the number of 
objects in each redshift bin. Error bars arc the scatter of the distributions. 
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the local Universe. From a Kuiper test, the field and cluster 
ETG population were shown to be different. 

We extend Kuiper and Kolmogorov Smirnoff tests anal- 
ysis to our group and field passive early-type galaxy size 
distribution in different mass bins to test environmental de- 
pendences on their distribution scatter that might distin- 
guish size evolution in the groups and the field. We divide 
our sample in 3 redshift bins (0.2 < z < 0.5, 0.5 < z < 0.8 
and z > 0.8) and 2 mass bins (10.5 < log(M/M Q ) < 11 and 
log(M/M @ ) > 11). 

A point to take into account when normalizing with the 
local relation is the fact that groups and field galaxies in the 
local universe could follow different mass-size relations. We 
have checked this point by comparing the medi an mass-size 
relatio n of early-type galaxies in groups f rom thelYang et alj 
(120071 ) catalog to the full DR7 sample. I Yang et alj (|2007l ) 
put together a catalog of ~ 300000 groups detected in the 
SDSS DR4 using an automated halo-based group finder. 
For this work, we restricted to groups with more than 2 
members and removed those objects affected by edge effects 
(fedge < 0.6). According to their figure 2, ~ 80% of the 
groups have ~ 20% contamination which is comparable to 
the expected contamination in our sample. We use as halo 
mass estimate, HM1, which is based on the characteristic 
luminosity of the group but results remain unchanged when 
using an halo mass estimate based on the characteristic stel- 
lar mass. The group catalog has been correlated with the 
catalog of 2D sersic decompositions by Simard et al. (2011) 
to get a size estimate for all our galaxies (circularized effec- 
tive radius of the best fit) as well as with t he morphological 
catalog by iHuertas-Companv et alj (|201ll ) to select ETGs 
(Pearly > 0.8). The mass-size relations that we obtained for 
the field and group galaxies are fully consistent within la, 
we therefore use the same local relation to normalize group 
and field sizes. 

Results of the analysis are summarized in table and 
the different size distributions for several sample selections 
are shown in Fig. 1121 For all redshift and mass bins, both the 
KS and Kuiper test results almost always give a probability 
> 80% that field and groups size ratios are drawn from the 
same distribution. 

Our main result is therefore that the mass-size relation 
does not depend on environment for field and groups with 



halo masses M2ooc/Mq < 10 



One might argue that 



this lack of dependence with environment could be a conse- 
quence of field contaminations in the group sample, specially 
at large distances from the group center. We have checked 
though that the distributions remain the same when con- 
sidering only galaxies with l arge probabilities to be group 
members (Pmem > 0.8 from lGeorge et al1l201ll ) and close 
to the center (d < 0.5 x R200) where we expect a contami- 
nation of » 15%. (see appendix [B}. 

Fig. [TjJ] also shows the evolution of the size of field and 
group passive galaxies divided into the same morphological 
and stellar mass bins than in the previous section. Relations 
are more scattered because of the Poisson noise (especially 
for massive lenticulars) but we do recognize the same trends 
as for the mixed population. 

Recent works at similar redshifts hav e found that 
galaxy size evolution depends on environment. ICooper et all 
(20 Hal ) find indeed that denser environments at z < 1 
tend to be populated by larger galaxies. Differences with 



this work can come from the way environment is measured 
and/or how galaxy selection is performed. Cooper and col- 
laborators selection is indeed based exclusively on the Sersic 
index, which as shown in Fig.[5]leads to the inclusion of some 
star-forming galaxies which tend to increase the scatter of 
the field population. As a matter of fact, we find that KS 
and Kuiper tests applied to the size distributions of group 
and field galaxies selected on the basis of the Sersic index do 
show that the distributions have probabilities < 50% to be 
drawn from the same distributions. However, median sizes 
do not change significantly so this fact cannot fully explain 
the difference between the two works. Concerning the en- 
vironment measurement, we use here the DM halo mass as 
primary environment estimator whereas Cooper and collab- 
orators use the local density based on neighbors. 



3.3.1 BGGs 

Brightest group galaxies (BGGs) deserve a particular men- 
tion since they have been more studied because of their very 
particular position at the center of massive structures and 
their high surface brightness which make them easy to detect 
and analyze. Current works, though, do not agree about the 
BCG size evolution and how BCG sizes compare to those of 
field and satellite galaxies. 

Based on SDSS data. lBernardl(|2009l ) found that BCGs 
are larger than field and satellite galaxies at fixed stellar 
mass and that there is a steep evolution of their size from 
z ~ 0.3 to present. The author also argues that minor dry 
mergers are the most probable mechanis m to explain the 
build- up of these objects. Also in the SDSS. IWeinmann et al.l 
(|2009h did not find a significant difference between the sizes 
of centrals and s atellites in groups. 

At higher z, lAscaso et all (|201lD find a significant size 
evolution from z ~ 0.6 to present, but do not detect any 
evolution in their profile. They interpret this fact as a sig - 
nature of feedback instead of merging. (|Stott et alj 12011 ) 
studied a sample of high redshift BCGs and found a very 
mild evolution of the BCG size from z ~ 1. 

From the modeling point of view, IShankar etakl (|201lh 
showed that BCGs should evolve much faster than satellite 
galaxies. 

We use the BGGs definition as the most massive galax- 
ies within an NFW scale radius of the X-ray position, from 
the George et al. (2011) catalog. Fig. 1 141 shows the size evo- 
lution for passive BGGs and satellite group members with 
similar stellar mass (log(M /Mq > 11). We find that the two 
population evolve in a similar way within the error bars, e.g. 
we do not observe significant differences in the size evolu- 
tion of BCGs as compared to satellite group members with 
similar stellar mass. 



4 DISCUSSION 

To understand why not all populations evolve in the same 
way, we consider different scenarios of galaxy growth. There 
are two main, well distinct, physical mechanisms proposed 
so far in the literature to puff up massive bulges from high 
redshifts to the local Universe: galaxy mergers and mass loss 
via quasar and/or stellar wind. 
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Figure 12. Size ratio distributions of group (filled histograms) and field galaxies (empty histograms) for different sample selections and 
different redshift and stellar mass ranges. The red dashed vertical line indicates a size ratio of 1. Each column shows a different rcdshift 
bin (0.2 < z < 0.5, 0.5 < z < 0.8 and 0.8 < z < 1.1 from left to right). Each row shows a different morphology and mass selection as 
explained in the text. 
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0.99 
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53 


89 






PEs 


0.95 


0.99 


31 


51 






PSOs 


0.55 


0.91 


22 


38 



Table 2. Results of Kuiper and Kolmogorov-Smirnoff statistical tests applied to group and field galaxies size ratio distributions in 
different samples (see text for details). Pall stands for all passive galaxies, PEs are passive ellipticals and PSOs passive lenticulars. N group 
and N fi e id indicate the number of galaxies in the group and field samples respectively. We did not compute the statistical tests when 
the number of objects was below 10. 



4.1 Mergers 

Hierarchical models of structure formation envisage the 
growth in size of massive galaxies via a sequence of major 
and minor mergers. While (mainly gas-rich) major mergers 
are believed to happen at high redshifts and may be re- 
sponsible for forming the galaxy, minor dry merger happen 
on cosmological timescales and tend to impact mainly the 
external regions of the galaxy, thus increasing its size, but 
leaving its inner regions mostly unaltered. 

In this work we consider the predictions of several, 
representative, hierarchical galaxy evolution models that 
differ in terms of underlying techniq ues and phys i cal as - 
sumptions. We adop t, more prec i sely, 



iHopkins etlll (|2009l ). iGuo et all (|2010l ) and IShankar et all 



Bower et all J200fj). 



( 20 111 ). All model predictions have been computed for the 
range of stellar masses Q and galaxy type of interest to this 
paper. We mainly considered galaxies with B/T > 0.5 when 
comparing to early-type galaxies (though we checked that 
our results are practically unchanged when restricting to 
B/T > 0.7), and galaxies with 0.3 < B/T < 0.7 when com- 
paring to SOs. 

Both the iBower et all (|2006t ) and IGuo et alj (|2010h 
models follow the hierarchical growth of gal axies along the 
merg er trees of the Millennium simulation (|Springel et al.l 
2005J). Galaxy progenitors are initially disk-like and after a 
major merger the remnant is considered to be an elliptical 
(though disk regrowth can happen). Minor mergers instead 
tend to preserve the initial morphology of the most mas- 
sive progenitor but tend to increase the mass of the bulge 
and disk components via the aggregation of old stars and 
newly formed ones during the merger. Half-mass sizes are 
then updated at each merger event assuming energy conser- 
vation. Despite being built on the same dark matter simu- 

2 All predictions have been corrected to a common Chabrier IMF. 
We however expect differences in the IMF to have a minor impact 
on the rate of size evolution of massive spheroids. 



lation, the subhalo/galaxy merger rates of these two models 
differ due to the different corrections in dynamical friction 
timescales. In both models, bulges can also grow via disk in- 
stabilities. However, the implementation of the latter phys- 
ical ingredient sub stantially differs in the two models, with 
I Bower et al. I (|2006h assuming a much stro nger bulge g rowth 
via disk instabilities wit h respect to th e Guo et all (|2010t ) 
model (see discussion in lShankar et al]|2012l ). Most impor- 
tantly for size evolution, both models do not consider gas 
dissipation during gas-ric h major merger s , a fea ture that has 
instead been included bylShankar et al] (|201ll) by properly 
adapt ing the I Guo et al] (|2010l ) model. The Hopkins et al.l 



(2009) model follows the analytic mass accretion histories 
of haloes and at each time step initializes central galaxies 
and infalling satellites according to empirical correlations 
inspired by halo occupation techniques and high-redshift 
data. Equivalently to the models discussed above, at each 
merger the half-mass radius is updated following energy 
conservation argum ents, with also gas dissipation (see also 
iNipoti et al.1 12012 | for a more recent work adopting simi- 
lar techniques). Hopkins et alj |2009l ) have mainly focused 
on early-type galaxies (not lenticulars) with stellar masses 



above M star > 10 M 



©■ 



4-1.1 The Size Growth of Ellipticals compared to model 
predictions 

In Figs. 1151 we compare our results with predictions from 
Shankar et al. (2011), Bower et al. (2006) and Guo et al. 
(2010). For ellipticals, we focus here on the lowest and high- 
est stellar mass bins of fig 1111 where the differences are 
stronger (i.e. 10.5 < log(M/M Q ) < 11 and log(M/M Q ) > 
11.2). The first relevant feature arising from the comparison 
is that most of the theoretical predictions are in agreement 
with the (mild) size evolution of ETG galaxies with stellar 
masses M* < 10 11 (left panel of Fig. I15[) (specially taking 
into account the uncertainties due to morphological classifi- 
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Figure 13. Evolution of the size of galaxies from our sample as a function of redshift referred to the local relation from SDSS (see 
text for details). Filled circles are group galaxies and empty squares are from the field. Error bars are related to the scatter of the size 
distribution (see text) and numbers indicate the number of galaxies in each redshift bin in the group sample. 



cations). Merger models are therefore successful to predict 
the size evolution of these lower mass ellipticals. 

More interestingly, despite the significant variance in 
the input parameters and/or physics most merger mod- 
els seem to be unable to reproduce completely the fast 



drop in sizes for quiescent ellipticals with stellar masses 
Log(M,/Mo) > 11.2, especially at 2 > 0.5 (right panel 
of Fig. I15[) . This is in line, and further complements, the 
recent claims of a poss ible inefficiency of the puffing-up via 
mergers pointed out bv lShankar et al.l (|201ll ), ICimatti et al.l 
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Figure 14. Size evolution of BCGs (filled circles) compared to satellite group members with log(M / Mq) > 11 (empty circles). Numbers 
show the number of BCGs in each bin. 



l|2012l ). and iNipoti et all l|2012l ). Notice that in this stellar 
mass bin, the contribution of SO galaxies is almost negligible 
(see fig. [Tl)) , so a similar behavior is found when selecting all 
ETGs with log(M,/M Q ) > 11.2 and should be independent 
of the morphological selection. 

INipoti et alj (|2012l ) have recently shown the results of 
a merger model initialized via halo occupation techniques 
that, by neglecting dissipative processes and assuming only 
mergers with spheroids, maximizes the evolution in surface 
density. They conclude that minor and major mergers may 
not be sufficient to ex plain the observed size grow th of early 
type galaxies (see also lLopez-Saniuan et al1 [2012). An other 
work where the size gr owth since z ~ 1 is ex plained almost 
completely by mergers (jNewman et alj|2012h has to assume 
very short merger time scales (1 Gyr) and steep growth effi- 
ciency, optimistic with respect to observations (notice in ad- 
dition that the size evolution they measure is less steep than 
the one reported here though, be cause of the select ion used, 
as shown in sections |3. II and l33)) . IOser etall (|2012h have re- 
cently analyzed 40 cosmological re- simulations of individual 
massive galaxies with final M* > 6.3 x 1O 1O M0, out of which 
25 appear quiescent early-type galaxies. While they claim a 
strong size evolution in the cumulative distribution of galax- 
ies with present stellar mass M» > 6.3 x 1O 1O M0, a close 
inspection of their Figure 1 (left panel) reveals that galaxies 
above M* > 2 x IO^Mq have indeed had a rather mild size 
evolution at fixed stellar mass of about ~ 30 — 40% at z < 1. 
Nevertheless, this may not necessarily reflect a failure of hi- 
erarchical models, as part of the discrepancy could simply 
arise from the specific underlying assum ptions made. For 
examp le, the early semi-analytic model bv lKhochfar fc Silkl 
(|2006h predicts a stronger size evolution for the very massive 



galaxies with M* > 5 x 10 11 xMq, in better agreement with 
observations. 

Overall, most merger models have some difficulties in 
fully reproducing the size evolution of the most massive 
early-type galaxies. Part of the apparent evolution in sizes 
may be driven by relatively younger larger galaxies formed 
at later epochs (progenitor bias) even if the relevance of 
this effect is still u n clear. A recent empirical model by 
lLopez-Saniuan et al.l (|2012l ). based on merger observations 
in the COSMOS field, shows in fact that taking into consid- 
eration major and minor mergers observed in COSMOS can 
explain ~ 55% of the size evolution of massive (> 10 Mq) 
ETGs. If the progenitor bias of massive ETGs accounts for 
a factor 1.25, this work can explain ~ 75% of the size evo- 
lution. We also show, for completeness, in figure [TS] the ex- 
pected evolution taking into account this effect from that 
work (notice however that the selection in terms of morphol- 
ogy and star-formation is not exactly the same than the one 
used in this work). However. I Whitaker et all (|2012l ) showed 
that the recently quenched galaxies at 1 < z < 2 are not 
significantly larger than their older counterparts, suggesting 
that the effect of the so called progenitor bias is limited. 
Moreover, fig. [9] also shows that the number densities (frac- 
tions) of massive ellipticals do not evolve significantly from 
z ~ 1, pointing again towards a reduced effect of newly 
formed galaxies. 

There are also other possible tensions between merger 
models and current observations. A puffing-up of the exter- 
nal regions wi th no inner densit y variation as suggested by 
mergers (e.g.. iNaab et al1l2009l ). naturally produces an in- 
crease in the Sersic index n of an ideally light profile fitted 
to the projected density profile. In Fig. 1161 we plot the me- 
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Figure 15. Size evolution of elliptical passive galaxies in two mass bins. Left panel: 10.7 < log(M / Mq) < 11, right panel: log(M / Mq) > 
11.2. Numbers indicate the number of objects in each redshift bin. Error bars are the scatter of the distributions. The different lines (see 
legend) show the prediction from different semi-analytic and empirical models (see text for details). 
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Figure 16. Median Scrsic index evolution for elliptical galaxies in the mass range 11.2 < log(M /Mq) < 12. The local point comes from 
the Scrsic decompositions of Simard ct al. (2011) on visually classified ellipticals by Nair & Abraham (2010). 
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dian Sersic index for massive elliptical galaxies. The median 
Sersic index of the overall population does not evolve since 
z ~ 1. These results are in possible tensi on with this predic - 
tion (similar conclusions were found by IStott et al.l (|201lf ) 
for instance). 

Another prediction of merger models is that galaxies 
residing at the centre of more massive haloes should, at 
fixed stellar mass, experience more mergers and thus be 
larger than their counterparts in less massive haloes , at 
least above M u > 5 x 10 12 Mq, according to the analysis of 
IShankar etafl (|201lh . The latter, in fact, showed that BCGs 
should evolve much faster at z < 2, and also end up being 
larger than other galaxies of similar mass. We do not observe 
this trend. 

4-1.2 Environmental dependencies compared to model 
predictions 

Concerning the dependence of galaxy sizes and their evolu- 
tion on environment, we compare our observations to pre- 
dictions from Shankar et al. (2011) in figure [T71 We plot the 
size-halo mass relation for two redshift bins from COSMOS 
(0.5 < z < 0.8 and 0.8 < z < 1.0) as well as for the local 
Universe from SDSS. 

Uncertainties in the observations have been calculated 
by bootstrapping 1000 times and recomputing the median of 
the distribution each time, as in all plots shown previously. 
To properly compare observations with model predictions at 
each redshift and halo mass interval of interest, we perform 
1000 Monte Carlo realiz ations in which we d raw subsam- 
ples of galaxies from the lShankar et all (|201 lh catalog with 
numbers equal to those in the SDSS and COSMOS sam- 
ples. Galaxies are selected to have B/T > 0.5 and to share 
the same stellar and halo mass intervals as in the observa- 
tions, a nd sizes have b e en no rmalized to the local mass-size 
relation IShankar et al.l (|201ll ). To each mock subsample we 
subst itute 30% (when comparing to COS MOS.lGeorge et all 
l201ll ) and 20% (when comparing to SDSS. lYang et alj|2007l ) 
of members with galaxies of the same stellar mass residing in 
the halo bin with the lowest mass in order to mimic contam- 
ination from the field. The latter has in general, as expected, 
an effect of partly reducing the increase of mean size with 
halo mass by up to 20%. For each Monte Carlo realization 
we compute the mean, and then extract, from the full dis- 
tribution of means, the final mean and its 1-sigma error. We 
do not consider errors in halo mass in our Monte Carlo real- 
izations. We believe that including them should increase the 
scatter among neighboring bins in halo mass and thus con- 
tribute to reducing the signal of increasing size with halo 
mass. However, we find from preliminary simulations that 
up to a factor of 2 statistical error in halo mass does not 
severely impact the final results. We will develop this point 
in detail in future work (Shankar et al. in prep.) 

As discussed in section 3 and summarized in Fig. 1131 
this work does not show any significant dependence of me- 
dian galaxy size on large-scale environment, e.g., galaxy sizes 
have similar medians in the field and in the groups, as well 
as in the SDSS groups from Yang et al. ( 2007) and the 
field. In contrast, Figure [17] shows that the I Shankar et al.l 
l|201ll ) model predicts a clear dependence of galaxy size on 
the host halo, predicting sizes in the SDSS and COSMOS 
group mass ranges to be around 1.5 times larger than those 



in the field. For the SDSS sample, this result is solid at 
> 3<r. The COSMOS sample shows the same trend, even 
if, because of smaller statistics, the disagreement between 
predictions and observations is only at ?a 1.5 — 2a. 

4-1.3 The Size Growth of Lenticulars compared to model 
predictions 

For what concerns lenticular galaxies fFig. H8|) . model galax- 
ies have been selected to have 0.4 < B/T < 0.7. For this 
range of B/T, disk instabilities also play a non-negli gible 
role in building bulges (see, e.g., IShankar et al.ll201ll . and 
references therein). Despite the significantly different proce- 
dures to grow bu lges in the models , from mild instab ilities 
|Guo et all |2010| ) to violent ones (|Bower et al.ll200rj ). the 
result in the median size evolution is similar to the ellip- 
tical results. It is interesting to note that, analogously to 
what inferred with respect to more bulge dominated galax- 
ies, hierarchical models are able to reproduce the trend in 
size evolution for lenticulars below M» < 10 11 Mq, but they 
tend to predict a shallower evolution above this mass. 

4.2 Expansion 

An alternative model for efficiently puffing up sizes of mas- 
sive, early-type galaxies, considers the galaxy expansion con- 
sequent ^_to_sjgjrificant nmss lo^ stellar 
winds ijFan et al] boost l2010l : iDamianov et alJl2009l ). This 
model envisages that all massive (M» > 3 x 10 10 Mq) early- 
type galaxies forming at z > 1 went through a rapid expul- 
sion of large amounts of mass, possibly caused by a powerful 
quasar wind, that caused the galaxy to expand. Numerical 
experimen ts in favor of this physical scenario h ave been per- 
formed bv lRagone-Figueroa fc Granatd (|201ll ). who showed 
that even in the presence of large amounts of dark matter, 
massive galaxies can significantly expand their sizes by a 
facto r approximately proportional to the fraction of mass 
loss. iFan et all (|2008h predict that most of the size evolu- 
tion for the massive spheroidal galaxies should be delayed 
with respect to the peak of quasar activity by about 0.5-1 
Gyr. The model also predicts a milder and possibly longer 
size evolution for early-type galaxies with stellar mass below 
M» < 2 x 10 10 Mq, for which the dominant energy input in 
the interstellar medium comes from supernovae explosions. 

Clearly, the typical evolutionary timescales to fully 
evolve a g alaxy onto the local size-mass relation in the 
IFan et al] (|2010h model are in general quite shorter than 
the cosmological ones required by a merger scenario. While 
a fast evolution in the size growth of massive ellipticals is 
supported by our data, their model predicts strong size evo- 
lution mainly at z > 1, while our data show strong size 
evol ution also at z < 1, at least for massive ellipticals. 

IFan et all (|2010h also point out that a fast size evolu- 
tion necessarily should lead to the co-existence of large and 
compact quiescent galaxies at any redshift z > 1, and thus 
a large dispersion up to a factor of < 5 — 6 in the size dis- 
tribution at fixed stellar mass. The dispersion should then 
significantly reduce below z < 1, as most of the galaxies 
should be already evolved and their formation rate, which 
paral lels the one of quasars peaking at z = 2 (|Lapi et al.l 
bOQrj ). should progressively drop at late times (see their Fig. 
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Figure 17. Size of elliptical galaxies as a function of halo mass for different redshift bins. Red diamonds are objects with 0.8 < z < 1.0, 
orange squares are 0.5 < z < 0.8. The black dashed-dotted line, the orange dashed line and the red dashed-thrce dotted line show the 
expected relation from Shankar et al. (2011) models at z = 0, z = 0.5 and 2 = 1. Values from models have been normalized so that, by 
definition, the field observed value at an halo mass of log(Mt l /M l r l ) = 12.8 — 13 is equal to 1. The solid black line and the filled squares 
are the values measured from the SDSS and the group catalog from Yang ct al. (2007). Errors in models and observations are errors on 
the median values computed through bootstrapping (see text for details). 



1). The dispersion in sizes for quiescent, massive early-type 
galaxies that we measure from COSMOS seems instead to 
be rather contained, within a factor of two at fixed redshift 
and stellar mass. We also note, however, that our lower mass, 
quiescent ETG that have larger sizes might be, at least in 
part, be composed by the population of high-r edshift, al- 
ready evolved early-type galaxies predicted by lLapi et all 
|2006l ). 

The expansion model predicts that compact galax- 
ies should be relatively young at the time of observa- 
tion, being close to their formation epoch (because ex- 
pansion happens less than a Gyr after the quasar epoch). 
This may be possibly at variance with the rather old 
ages that usually characterize mas sive ellipticals (see also 
iRagone-Figueroa fc Granatd 1201 ll for similar considera- 
tions). 

On the other hand. lRagone-Figueroa fc Granatol (|201ll ) 
also found in their numerical tests that the galaxy mass 
profiles should not change after the blow-out, if the mass 
loss is contained to a factor of two or so. This might explain 
the nearly constancy with redshift of the Sersic index in 
our sample. Also, being galaxy expansion an in-situ physical 
mechanism, it should be largely independent of environment, 
as suggested by our data (Figs. [131 and [17)). 

In summary, while merger models can explain most of 
the size growth size z ~ 1, we find evidence that both pro- 
posed scenarios (mergers and mass loss via quasar and/or 
stellar winds) suffer from some shortcomings with respect 



to our results and more advances in modeling are clearly 
needed to deeply understand how the two different scenario 
shape the galaxy size evolution. 



5 SUMMARY AND CONCLUSIONS 

We have studied a sample of 3146 group galaxies and 3760 
field galaxies from the COSMOS survey, and selected 298 
group and 384 field galaxies as passive galaxies (based on 
their M(NUV)-M(R) dust corrected rest-frame color) with 
log(M /Mq) > 10.5. We show, for the first time, how the 
mass-size relation and size growth depend on the detailed 
morphology of the quiescent population (mainly ellipticals 
and lenticulars) as well as on large-scale environment defined 
by the dark matter halo mass. 
Our main results are: 

(i) A detailed morphological dissection of the passive pop- 
ulation up to 2 ~ 1 reveals that: 

• About 80% of all passive galaxies have an early-type 
morphology at all stellar masses and at all redshifts from 
z ~ 1. The remaining 20% are essentially early- type spi- 
rals. 

• Early-type galaxies are both ellipticals and SOs. At 
all redshifts up to z ~ 1, elliptical galaxies dominate 
the ETG population at masses log(M/M Q ) > 11 - 11.2. 
At z < 0.5, galaxies with masses log(M / 'Mq) < 11 are 
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Figure 18. Size evolution of lenticular passive galaxies in two mass bins. Left panel: 10.5 < log(M / 'Mq) < 11, right panel: log(M / Mq) > 
11.0. Numbers indicate the number of objects in each redshift bin. Error bars are the scatter of the distributions. The red dashed line in 
the left panel indicates the median evolution of star-forming disks in our sample divided by a factor of 1.3 to reach the same normalization. 
The other lines (see legend) show the prediction from different semi-analytic models (see text for details). 



around half ellipticals and half lenticular, while lenticulars 
dominate the low-mass fractions at z > 0.5. An impor- 
tant fraction of group and field low-mass passive galaxies 
in the redshift range 0.5 < z < 1 are lenticular galaxies, 
e.g. have a disk component (see also Bundy et al. 2010, 
iMei et al.ll2012i ) 

Therefore, studying the population of passive galaxies as 
a whole mixes different morphological populations which 
do not necessarily share the same evolutions. 

(ii) When separating the ellipticals from the lenticulars, 
we show that galaxy size evolution strongly depends on mass 
and ETG morphology: 

• Massive ellipticals (log(M, /Mq) > 11.2) do experi- 
ence a very strong size evolution from z ~ 1 to present 
(r e oc (1 + z )-°-9S±o-i8y Even though the trend is some- 
how steeper than the one predicted by most of published 
semi-analytical models, minor dry mergers remain the 
most plausible explanation to the expansion. However, the 
Sersic index of this population is not significantly evolv- 
ing with time as expected from hierarchical models. We 
cannot exclude from our data that expansion via feed- 
back might also play a role on the size evolution of this 
population. 

• Less massive ellipticals (10.5 < log(M* /Mq) < 11) 
do not evolve much from z ~ 1 (from ~ 10% to ~ 30% 
depending on which morphological classification we use). 
This behavior is well reproduced by most of current semi- 



analytic models from which the size growth is mainly due 
to major and minor mergers. 

• The evolution observed in the lenticular population 
does not change significantly with stellar mass and they 
do show a size growth of 55% ± 10% since z = 1. 

(iii) Finally, we studied environmental effects on the size 
evolution by dividing our sample into field and group galax- 
ies: 

• We do not detect any significant evidence that 
the evolution of field and group galaxies (13 < 
log(M h /M ) < 14.2) with stellar masses log(M t /M e ) > 
11.2 are different. This is observed both i n the local 
Unive rse, when comparing SDSS groups from I Yang et al.l 
(I2007T I with field galaxies, and in our COSMOS group and 
field samples. This is also true for massive central ellip- 
ticals, where models cl e arly p redict a major impact from 
mergers IShankar et al.l |201ll ) . 

This result is at va riance with the stand ard hierarchi- 
cal model (e.g., from IShankar" et alJ (|201lh ), that pre- 
dicts that galaxy sizes in our group mass range should 
be « 1.5 larger than in the field at all redshifts. For the 
SDSS sample, the difference between median sizes from 
the observations and the model predictions is larger > 3<r. 
For the COSMOS sample, we observe the same trend at 
« 1.5 — 2<r, because of smaller statistics. 

• BGGs and satellite galaxies of similar stellar mass 
(log(M*/Mo) > 11) evolve in a similar way. This is in 
disagreement with hierarchical models, that also predict 
a difference between centrals and satellites living in ha- 
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los with similar masses than our groups (log(M, /Mq) ~ 
13.4— 13.6). In fact, Shankar et al. (2011) has shown that 
galaxies residing at the centre of more massive haloes 
should, at fixed stellar mass, experience more mergers 
and thus be larger than their counterparts in less mas- 
sive haloes. 

Future work will involve a detailed detailed discussion 
of environmental effects on models (Shankar et al. 2012, in 
prep) as well as an extension to higher redshifts by studying 
the sizes of ETGs in distant clusters of galaxies (Delaye et 
al. 2012, in prep). 
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APPENDIX A: VISUAL MORPHOLOGICAL 
CLASSIFICATIONS 

Some results presented in this paper are based on an auto- 
mated morphological classification of galaxies at high red- 
shift. Despite of the extensive tests performed to properly 
assess the accuracy of our classification, the separation of 
ellipticals from lenticulars might still remain a challenge. In 
this appendix, we want to make sure that our results do not 
change while using an alternative classical visual classifica- 
tion. Two of us (SM and MHC) performed independent de- 
tailed visual classifications of all passive galaxies previously 
selected as early-type by galSVM (we assume that the sepa- 
ration between early and late is good enough). SM classified 
galaxies into three morphological classes (Ellipticals, Lentic- 
ulars and early-spirals) while MHC followed a slightly differ- 
ent criterion and separated galaxies into two classes (disk, 
no disk). In Figs. IA2l and lA3l we reproduce figures ITUl and [TT1 
of the main text using the 3 morphological classifications (2 
visual and 1 automated). While there are some differences 
between the three estimates, the main trends discussed in 
the text do hold and the results are consistent within the er- 
ror bars. This confirms that our morphological classification 
is robust and that the major results presented in the paper 
are not biased by the automated classification. Most of the 
differences between visual and automated classifications are 
seen in the low mass elliptical population for which visual 
classifications show a steeper evolution than that estimated 
while using an automated classification. When the evolution 
of low mass ellipticals is compared to their massive counter- 
part, it is still less steep as discussed in section 13.2.11 and 
shown in table |A"T1 



APPENDIX B: ENVIRONMENT AND GALAXY 
SELECTION 

We include here the equivalent of Fig. [T^] but only with 
galaxies having a probability Pmem > 0.8 to be a group 
member and being at a distance d < 0.5 x -R200 of the cluster 
center. Our results do not change and remain basically the 
same as the ones shown in table [2] The size distribution in 
mass and redshift bins does not depend on environment. 
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Figure Al. Morphological fractions of passive galaxies as a function of stellar mass and redshift using the visual morphological clas- 
sifications by MHC (solid lines) and SM (dotted-dashed lines). For MHC: red solid lines are ellipticals and blue are galaxies with an 
observed disk component. For SM: green lines are Sas, cyan lines are SOs and red lines are ellipticals. The dashed-dottcd blue line shows 
the fraction of galaxies with an observed disk component (SOs+Sas). The uncertainties are calculated following Gehrels (1986; see Section 
3 for binomial statistics; see also Mei et al. 2009). These approximations apply even when ratios of different events arc calculated from 
small numbers, and yield the lower and upper limits of a binomial distribution within the 84% confidence limit, corresponding to la. 



Stellar mass bin 


Morphology 


a 


10.5 < log{M/M Q ) < 11 


galSVM 


0.34 ±0.17 




MHC 


0.52 ± 0.23 




SM 


0.59 ± 0.27 


11 < log(M/M Q ) < 11.2 


galSVM 


0.63 ± 0.21 




MHC 


0.85 ± 0.22 




SM 


0.61 ± 0.23 


11.2 < log(M/M e ) < 12 


galSVM 


0.98 ± 0.21 




MHC 


0.94 ±0.21 




SM 


0.95 ± 0.21 



Table Al. a values parameterizing the size evolution of ellipticals (r e oc (1 + z) a ) in three mass bins obtained using three different 
morphological classifications. 
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Figure A2. Size evolution of passive ellipticals (top row) and SOs (bottom row) in two different mass bins using three different mor- 
phological classifications (see text for details). Empty circles show the automated classification from galSVM, squares are the visual 
classification from SM and triangles indicate are obtained with the visual classification of MHC. 
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Figure A3. Size evolution of elliptical (left column) passive galaxies in three mass bins using three different morphological classifications 
(see text for details). Empty circles show the automated classification from galSVM, squares are the visual classification from SM 
and triangles indicate are obtained with the visual classification of MHC. Top row: 10.5 < log(M /Mq) < 11, middle row: 11.0 < 
log(M I 'Mq) > 11.2 and bottom row: 11.2 < log(M j 'Mq) > 12. Numbers indicate the number of objects in each redshift bin. Error bars 
are the scatter of the distributions. 
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Figure Bl. Same as figure [T2l but only for galaxies having a probability Pmem > 0-8 to be a group member and being at a distance 
d < 0.5 X -R200 of the cluster center. Wherever only one histogram is shown it is because there are too few objects in that bin. 



